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Dynamical Virial Relations 
and Invalidity of the Boltzmann Kinetic Equation 

Yu. E. Kuzovle\0 
Donetsk Institute for Physics and Technology of NASU, 83114 Donetsk, Ukraine 

A sequence of exact relations is found which connect one- and many-particle time-dependent 
distribution functions of low-density gas with their derivatives in respect to mean density. It is 
shown that, at least in the context of spatially non-uniform gas evolutions, these relations forbid the 
"molecular chaos propagation" and imply inapplicability of the Boltzmann kinetic equation even 
under the Boltzmann-Grad limit and regardless of degree of the non-uniformity. 
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1. Introduction. 

Subject of our interest here will be physical status of 
the celebrated L. Boltzmann's kinetic equation (BE) for 
classical gas [l|-|ff . It can be written as 



F(t, r, v) = -vVF(t, r, v) + 
+ n C 2 F(t, r, v) * F(t, r, w) 



(1) 



where F(t, r, v) is one-particle distribution function nor- 
malized to volume, n is mean gas density, and C 2 is 
the Boltzmann's pair collision operator ("collision inte- 
gral" ) which acts onto F 's velocity argument as 



C 2 F(v) * F(w) = J d 3 w \v-w\J d 2 b x 
x [F(v in (b,v,w))F(w in (b,v,w)) - F{v)F{w)} 



(2) 



with b being the impact parameter vector (perpendicu- 
lar to v — w ) and v m and w m input velocities what 
lead to the given output ones. 

Anybody can agree that today's kinetic theory, as well 
as 100 years ago, is unthinkable without BE. Neverthe- 
less, nobody have presented a rigorous and at the same 
time general enough derivation of BE from such exact 
equations of statistical mechanics as the BBGKY equa- 
tions (BBGKYE) 0^. 

In general, BBGKYE can only produce something like 



F = -vVF + nC 2 F*F + n 2 C 3 F*F*F- 
+ n 3 C 4 F * F * F * F + . . . 



(3) 



(see below). Therefore, BE follows from BBGKYE un- 
der formal "low-density limit" only, when n — > 04al • 
But this is non-physical limit, since it enforces the mean 
free path of gas atoms, A ~ l/i:a 2 n (with a denoting 
characteristic radius of atom- atom interaction), to tend 
to infinity. 

Much more physically reasonable idealization is the 
"Bpltzmann-Grad limit" (BGL) when n — >■ 00 and si- 
multaneously a — > , in such way that A is kept con- 
stant. At that, a 3 n — > , i.e. gas becomes "infinitely 
rare" . This gave rise to hopes that the terms with C3 , 
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C4 , etc., in EqJ3]must vanish under BGL too. Moreover, 
O. Lanford and others even suggested a proof of this hy- 
pothesis for hard-sphere gas 0, Q. But their interpre- 
tation of BBGKYE in this singular case was definitely 
wrong (in fact they from the very beginning substituted 
interaction terms there by pair collision operators) [9( . 

In any case, contributions with C3 , C4 , etc., in Eqj3] 
appear due to violation of the Boltzmann's "molecular 
chaos hypothesis", that is due to statistical correlations 
between two atoms what are in mutually input (pre- 
collision) states (approaching each other). Naturally, 
such correlations take place mainly at inter-atom separa- 
tions comparable with a , but this does not mean that 
they vanish under BGL (merely they sit just where they 
are most "harmful" ) . 

What is physical meaning of the pre-collision correla- 
tions ? 

Notice that BE (HJ-© fully corresponds to the 
"probability-theoretical view" at physical world: it pre- 
sumes that gas evolution consists of elementary random 
events (collisions with various input velocities and impact 
parameters) which occur independently one on another 
and all possess strictly certain (though may be unknown 
numerically) a priori (conditional) probabilities. 

Such the view was originated by J. Bernoully 300 years 
ago [l3| • But in 20-th century N. Krylov [llj showed that 
assumption that any sort of events can be furnished with 
certain a priori probability, or relative frequency, has no 
support in rigorous statistical mechanics. In opposite, 
mechanics allows different relative frequencies on differ- 
ent phase trajectories of a many-particle system. Then, 
fluctuations of relative frequencies from one phase trajec- 
tory to another, - being averaged over statistical ensem- 
ble of trajectories (experiments), - produce correlations 
between events and particles even in absence of cause- 
and-consequence connections between them. 

Just such statistical correlations were found and inves- 
tigated in my works [12| (see also [13|]) and [lJ-OJj on 
spatially inhomogeneous evolutions of fluids and random 
walks of their particles. And just such correlations de- 
termine higher-order terms in Eq|31 impeding reduction 
of BBGKYE to BE even in spite of BGL. 

At that, role of the inhomogeneity is to visualize and 
magnify effects of the uncertainty of collision probabili- 



ties of fluid's particles. In particular experiments it mani- F(t, r, v) and many-particle distribution functions (DF), 

fests itself in the form of scaless (1/f-type) low- frequency we will derive a sequence of exact relations between them 

fluctuations of both individual particles' mobilities and and their derivatives in respect to density n , and then, 

collective transport characteristics of a fluid. with the help of these relations, demonstrate that the 

tt T , r r • i-j-j. r i-u t^t^ prc-collision inter-particle statistical correlations stay fi- 

Here, i suggest new proofs of invalidity ol the BF *\ . .„ , „„ T 

0, . nl . , 11 TXTl c , ,. mtc and sigmflcant even under BGL. 

in spatially mhomogeneous problems [20|. starting b 

from general properties and exact density expansions of 
I 

2. Many-particle statistical dynamics . 

Let N 3> f gas atoms are contained, by means of an auxiliary external potential, in volume $7 = N/n . Let 
Xi = {ri,Vi\ denote variables of «-th, atom, 

Fi(t) = F(t,Xi) = Fi(t,Xi) , Fij(t) = F 2 (t,Xi,Xj) , F ijk (t) = F 3 (t,Xi,Xj,Xk) , ■■■ Fi,.. N (t) = 

Ff{(t, Xi,X2, . . . Xn) are one-, two-, three-, ... TV -particle DFs, 

Li , Lij , Lijk , ... Li...n are one-, two-, three-, ... N -particle Liouville operators (including the auxiliary 
potential), and 

S,(i) = exp(Li t) , Sij(t) — exp(Lijt) , ... Si...jv(t) = exp(£i...jvi) are corresponding evolution operators. 

All these objects are completely symmetric functions of their arguments (indices). 

If all DF are thought normalized to volume and all obey the requirement of mutual consistency, then 

1 '"^....(t) = 1 , Fi ....(*)= k / F 1 ... ss+1 (t) = -^- [ ... [ F lmN (t) , ( 4 ) 



" J\ J a i: Js+1 " Js+1 J N 

where J, . . . = L, d 3 r^ J d 3 ^ .... Evolution of all marginal DFs is determined by that of the whole system's DF 
according to 

Fi... N {t) = S 1 ... N (t)F 1 ... N {0) (5) 

Since it conserves full phase volume, the equalities ((H) will be satisfied at all times if it is so at one, "initial" , time 
moment, e.g. at t = . But some specific form of DFs can realize at one moment only. For example, we may choose 

*i....(0) =nj=i^(0) , (6) 

thus assuming "molecular chaos" at initial time moment. 

Further, it is convenient to introduce operation o of "coherent product" of the evolution operators. By definition, 
for any two non-intersecting sets of indices, 

Si...j o Sk...i — Si...jk...i (7) 

In essence, the left side here is mere equivalent notation for the right-hand side. For intersecting sets, of course, one 
must take their union. Thus, for instance, Si...j ■ — Si o . . . o Sj , and we can write identities 

N N 

si...n = Si.... n ° s i = Si - s n °t i + ^j- 1 )} = ( g ) 

j = S+l j = 8+1 

N 

= 5i.„ a + Y, S i-- ° ( S i ~ !) + 12 S i-s ° ( S i ~ !) ° ( Sk - !) + 

.7 = 8+1 s+l<j<fc<JV 

AT 

+ Yl 5 i-« ° ^ - 1) o (5 fc - 1) o (ft - 1) + . . . + Sl.., J] o(5, - I ) = 

s+1 <j<fc<( <JV i = s+l 

TV 

= Si. ..8 + 2_^ [Si. ..SJ — Si. ..8 ] + 2_^ [Sl...sj k — Sl... s j — Sl...sk + Si. ..s ] + ■ ■ • , 

3 = s+l s+l<j<fc<JV 

where < s < N (at s = , of course, both Si... s and Fi... s should be replaced by I). 
Combining these identities with Eqs|4]and[5] after simple algebra one obtains 



Fl .. A t) = si... s{ t)F,.,(o) + e kl J N :% k J s+1 ---J s j s ^m °[w*)-i]}^ 



Hfc(O) (9) 



fc = 1 ■-" >."■ "- -8fl «*+* j=1 



Advantage of this representation of the marginal DFs' evolution, in comparison with that given by (|4|) and (0, is 
that it can be directly extended to the thermodynamic limit, when N — > oo and ft — » oo (the auxiliary potential 
removes) at fixed N/il — > n = const . 

3. Density expansion. 

The Eq|S]by itself does prompt conditions of this extension. Namely, a limit form of the initial DFs -Fi... s (0) at 
£1 — > oo should be such that all the integrals in Eq[9]turn to well converging "proper integrals". That is such that, 
at any s and k , 

[S s+ i(t) - l]o ... o [S s+ k(t) - 1] i f i... ss +i... s +fc ->■ when r s+ i . . . r s+k -> oo , (10) 

in fast enough (integrable) way. This means that "at infinity" (asymptotically) all the initial DFs are invariant in 
respect to both translations and collisions of atoms. It will be the case if our gas is thermodynamically equilibrium 
(hence, homogeneous and rest) at infinity, i.e. 

Fi „.._i a -> Fi..,-iF (e «»(i; s ) at r s -> oo , (11) 

where F^ eq \v) = (27rT/m) _3//2 exp(— mv 2 /2T) is the Maxwell velocity distribution with some temperature T ( 
m is atomic mass). This agrees with requirements (j4]) in their limit form, that is, - in the case ([6]), - with 

lim - I d 3 n I (fiviFi = 1 (12) 

£Woo il J n J 

We see that the thermodynamic limit in Eq|9]presumes that our gas (or, generally, fluid) is essentially inhomogeneous 
(if not equilibrium at all) . 
Then EqJS] transforms to infinite series 



oo k p p 

Fi....(t) = 5i....(t)Fi....(0) + E TT / •/ { S i-At)H lS s+j (t) -I}} F u .. a+k (0) = 



(13) 



= Si... a {t) Fi... t (0) + n / [S 1 ... ss+1 (t)-S 1 ... s (t)}F 1 ... ss+1 {0) + 

Js+l 

"*" "of / / [ft-«8+ls+2(t) — 5*1 ...8S+l(t) — Si ...ss+2{t) + 5 1 ! ... s (t) ] F\ ...ss+ls+2(0) + ••• 

^ ! Js+l Js+2 

with initial DFs (and, hence, arbitrary-time DFs) satisfying conditions (fTQ|) and (fTTj) . and J = f d 3 rj J d 3 Vj .... 

The Eq |13l represents density expansion of total evolution operator of the infinitely-many-particle system. If its 
initial state is treated irrespective to n , then Eq 11 31 becomes density expansion of its time-dependent future (or past) 
DFs. This formal treatment is physically adequate at least when EqfT3] applies (as below) to sufficiently rare gas, in 
particular, characterized by initial "molecular chaos" ((6]). 

4. Dynamical virial relations . 

Below, we will confine ourselves by dilute gas, starting from the "molecular chaos" and being interested in the issue 
of its "propagation" with time. 

At s = 1 and s = 2 Eq[l3j as combined with (|6|), yields 

Fi(t) = Si(t)Fi(0) + nj[Si 2 (t)-S 1 (t)]Fi(0)F 2 (0) + (14) 

+ ^[[[Si 23 (t)-Si 2 (t)-Si 3 (t) + Si(t)]Fi(0)F 2 (0)F 3 (0) +... , 
F X2 {t) = S X2 (t)-Fi(0)F 2 (0) + n f[S 123 (t)-Si 2 (t)]F 1 (0)-F 2 (0)F 3 (0) + (15) 

+ ^ / /[ 5i23 4 (*)- S 123 (t) - S 12i {t) +S 12 (t)]F 1 (0) F 2 {0) F 3 {0) F 4 (0) + ... 

There is definite relation between these expressions. To perceive it, let us 
(i) introduce functions 



s oo fc ~ p s+k 

G 1 ... 8 {t) = \[o{s J {t)-i]F 1 ...M+Y,iA ••■/ n* 

3 = 1 k = l • Js + 1 Js+kj^l 



s+k 

,(t)-l]F 1 ... s+k (0) (16) 



(ii) return to EqfT3l expend there factors Si.., s (t) into "coherent product" of 1 + [Sj(t) — 1] 
(iii) take into account identity 



n k 



fc = i 



°= E IT/-/ iSi(t)-l}o...o[S k (t)-l}F 1 ... k (0) (17) 



what follows from Eqsl^land ll3l at s — (and seems trivial in view of general properties of Liouville operators), and 
(iv) with the help of (fTT|) rewrite F\ (t) and F\% (t) in the form 

F 1 (t) =F 1 (0) + G 1 (t) , n . 

F 12 (t) = F 12 (0) + Fi(0)G 2 (i) + F 2 (0)G 1 (t) + G l2 (t) {W> 

Next, let us apply here the initial "molecular chaos" (|6]) and consider the initial one-particle DF (normalized to 
volume according to p^j) 1 ) as independent on mean gas density n . Then, - treating (fTo| as functions of n too, - we 
obviously can write exact relations 

^ Gi... B (Jt) = J G u .. ss+1 (t) (19) 

Besides, differentiations of (JT7J) in respect to n produce identities 

J 1 ...J s G 1 ... s (t)=0 , / 1 .../ fc {n;=io[5 i («)-l]}nf=ii i i(0) -0 (20) 

At s = 1 equality (fT9|) . - together with EqsfT8l and [27)1 - gives 

^ Fi(t) = f[F 12 (t) - F 1 (0)G 2 (t) - G 1 (t)F 2 (0) - F 1 (0)F 2 (0)] = 
= f[F 12 (t)+G 1 (t)G 2 (t)--F 1 (t)F 2 (t)} = f [F l2 (t) - Fi(t)F 2 (t)} (21) 



Formulae ()T9|) and (|2Tj) are direct analogues of the "virial relations" found and considered in [15Hl9j | as inportant 
statistical properties of "molecular Brownian motion" . Here, I added the word "dynamical" in order to underline their 
principal difference from the relations known in equilibrium statistical mechanics. Another form of the "dynamical 
virial relations" (DVR) (fT9")l is expounded in Appendix 1 below. 

5. Symbolic kinetic equation . 

Kinetic equation (KE) is an equation what closely expresses, - like Eqs[TJand[3l - time derivative of the one-particle 
DF through it itself [2|. To construct a KE, first let us write 

L\... s = ii... s + L[ 



. S 7 



rO _ v^ rO Ti _ v^ ri (22) 

^l...s — 2^il<i<s n i ' ■'-'I... a ~ Z-*il<i<j<s n ij ' 

where operators L° and and V represent free motion of atoms and their (pair) interactions, respectively: 

L\ = - Vi V, , L' l3 = m- 1 *'( ri - r 3 ) {d/dv. L ~ d/dvj) , 

with $(p) being (short-range repulsive) interaction potential and $'(p) = d<fr(p)/dp . At once, introduce free 
evolution operators and "scattering operators" : 

S?.„,(t) = exp(L?... s t) , Z x ... s {t) = Si.„.(t)5?...,(-t) (23) 

(taking in mind that for any of our evolution operators 5 _1 (i) = S(—t) ). 

Then, perform time differentiation of EqUJJ] Using the boundary conditions at infinity, (fTU]l or (fTTj) . and the 
symmetry of DFs, is not too hard to make sure that, naturally, the result is the BBGKYE: 



L 1 ... s F 1 ... s (t) +^> / L' js+X F 1 

.7 = 1 Js+1 



F 1 ... a (t) = L 1 ....F 1 ... a (t)+y j n I L' ia+1 F 1 ...^. 1 (t) (24) 

Hence, we have only to express F\ 2 in terms of F\ 



With this purpose, let us treat Eg fMl) as series expansion of nF\{t) over ri-Fi(O) and try to invert it: 

Fx(0) = S 1 (-t)F 1 (t) - nS 1 (-t)f[S 12 (t)-S 1 (t)]S 1 (-t)S 2 (-t)F 1 (t)F 2 (t) + ... (25) 

Inserting this into Eqll5[ after a tedious non-commutative algebra, we obtain an infinite series whose first two terms 
look as follow, 

oo 

F 12 (t) = £n fe F«(*) = ( 26 ) 

fe = 

= Z 12 (t) Fi{t) F 2 (t) + nf[Z 123 (t) - Z 12 (t)Z 13 (t) - Z 12 {t)Z 23 (t) + Z 12 (t)}F 1 {t) F 2 (t) F 3 (t) + ... 

Here, a term with s-th degree of one-particle DF is of (s — 2)-th order in respect to density and contains s — 2 
integrations. Insertion of this expansion into Eq l24l at s — 1 yields the Eqj3]with 

C 2 F * F = J L' l2 Z 12 (t) F * F , (27) 

f [ L' 12 [ Z 123 (t) - Z 12 (t) Z 13 (t) - Z 12 (t) Z 23 [t) + Z 12 {t)] F * F * F (28) 

J 2 J 3 



C 3 F*F*F 



Spending a lot of time, one could continue series in (f2"5j) and (|2^|) and find also F^ 2 {t) and C4 . 

Notice that two written out right-hand terms of Eql26l confirm the results of semi- heuristic considerations [2. V&. l26j . 
although with essential differences. Namely, in [4| the second term is thought as a small correction to the first one, 
in the framework of the low-density limit (see Seel above), and time arguments of the "scattering operators" there 
are not unambiguously defined. In contrast to it, time argument of our operators Zij(t) , Z\ 23 (t) ,etc, is definitely 
total evolution time, so that they represent the whole pre-history of collisions. 

This difference is especially important from viewpoint of higher-order terms of Eg 1261 as considered in the framework 
of such more adequate approximation as the BGL. Actually, one can see that lower-order terms in Eg 1261 and hence in 
Eq|3] arise fully independently on higher-order terms in EqsQ~4l [T5l and [25l Therefore, if higher-order terms of Eol26l 
and [3] were insignificant then this would mean that arbitrary long evolution of one-particle DF is determined by only 
a few collisions. Absurdity of this enforced deduction clearly shows that in fact any of (infinitely many) terms of 
EqJ2S] essentially contributes to EqJ3]even under BGL. Hence, Eq(3]hardly can be useful in practice and sooner has a 
symbolic meaning only, while BE (TTJ) is invalid at all. 

Next, this logical necessity will be sustained mathematically. 

6. Violation of "molecular chaos propagation". 

The first term of the series EqUU practically reproduces the Bogolyubov's formulation [2J of the "molecular chaos" 
(MC) hypothesis: F\ 2 (t) => F^ 2 (t) — Z\ 2 {t) F\(t)F 2 {t) . Or, exploiting the pair correlation function, C\ 2 (t) = 
F 12 {t) - F l {t)F 2 {t) (see Appendix 1), C 12 {t) ^> [Z 12 {t) - 1] Fi(t)F 2 (t) What does it say? 

By definition of the scattering operators, 

Zi 2 (t)F(r x ,vi)F(r 2 ,v 2 ) = F(r' 1 + v' 1 t,v , 1 )F(r' 2 +v , 2 t,v' 2 ) , (29) 

where the primed variables represent such past state, time t ago, which lead to the given current state. Therefore, 
the MC hypothesis, - i.e. assumption about insignificance of second and higher terms of Eql^H - impliess that C\ 2 {t) 
differs from zero for two sorts of states. First, for currently interacting atoms, at \r 2 — r\\ < a . Second, for atoms 
in mutually post-collision (out-) states, at vi 2 ■ r\ 2 > , a < |ri2| < \vi 2 \t and | fo | < a. , where r% 2 = r 2 — J*i , 
v\2 = v 2 — vi , and b — r\ 2 — v\ 2 (y\ 2 ■ ri2)/|wi2| 2 is impact parameter vector (already mentioned in Seel). 

Thus, under the MC there are no correlations between atoms in mutually pre-collision states, - defined like the 
post-collision ones but with v\ 2 ■ r\ 2 < , - and no correlations between mere close though non-interacting atoms 
(for which |r"i2 1 is greater than a but comparable with a), but the payment for such pleasure is presence of the 
unreservedly far propagating post-collision correlations. 

If it was really so under BGL, then the exact relation (|2"Tj) . - after its multiplying by n , - would reduce to 

n^ Fi(t) = nf[Z l2 {t)-l]F 1 {t)F 2 {t) (30) 

The multiplication ensures finiteness of both sides here under transition to BGL, along with transition from n to 
physically more meaningful variable like k = ira 2 n = 1/A , so that nd/dn => nd/dn . 



In combination with Eqj29]this equality yields 

n— Fiftn.Vi) =n dr d 3 v 2 d 2 b\v 1 -v 2 \[F 1 (t,r 1 + {v' 1 -v 1 )T,v' 1 )F 1 (t,r 1 + {v' 2 ~v 1 )T,v' 2 ) 



-Fi(t,ri,vi)Fi(t, n + (v 2 -vi)t,v 2 )} = 

/ dr eiop(-TViV 1 )C 2 Fi(t,ri +vit,vi) Fi(t,ri +v 2 t,v 2 ) (31) 
Jo 



Here ' plays the same role as in the BE ([1}. And, as in ((T|), transition to BGL allows to neglect F\ 's changes at 
time and space scales related to a . 

Simultaneously, according to the previous section, the same MC assumption, - that Fi 2 (t) — Z\ 2 [t) F\{t)F 2 (t) , - 
produces the BE |T]) itself. 

Hence, if this is true assumption then Eqs l31l and [Tl should be compatible one with another. In fact, however, they 



can not be satisfied simultaneously, except purely spatially homogeneous case when Vi-Fi(i) = ! [27[ To become 
convinced of this fact in detail, one may e.g. consider linearized evolution of small local or periodic perturbations of 
equilibrium state ( F eq {v\) from above). 

Consequently, contribution of the rejected terms of Eq[26] into integral in DVR ([2~Tj) is on order of its value (see 
Eq |30[) under the MC assumption, and the latter, as applied to spatially inhomogeneous evolutions, is incompatible 
with the exact relation (|2"T|) even in spite of the BGL. 

In other words, the "molecular chaos propagation" fails in spatially inhomogeneous case, so that the BE proves to 
be invalid even under BGL [2a.|29j. 

I would like to underline that all the aforesaid equally embrace the hard-sphere interaction. 

To finish the paper and exclude hopes to "save" BE, we will supplement just presented proof of its invalidity with 
short notes on the pre-collision correlations. Additional comments on related many-particle correlations and their 
influence onto F\(t) 's evolution are placed to Appendix 2. 

7. Excess and pre-collision correlations. 

Let us return to Eq|26l and consider functions 

AF 12 (t) = Z 12 (t) AF[ 2 (t) = F 12 (t) - Z 12 (t) Frit) F 2 (t) = £~ =1 n k F$(t) = m) 

= nf 3 [Z 123 (t) - Z 12 (t) Z 13 (t) - Z 12 (t) Z 23 (t) + Z 12 (t) ] Fi (t) F 2 (t) F 3 (t) + . . . V ' 

They characterize those part of statistical correlations between two atoms what is due to not their interaction between 
themselves but common pre-history of their interactions (collisions) with the rest of gas. At that, k -th term of the 
sum represents connected chains (clusters) of at least k + 1 collisions (actual or virtual ones) conjointly involving 
k + 2 atoms. Statistical meaning of such "excess" (or "historical" [13]) correlations was exhaustively explained in 
QJ, EH (see also Seel above and notes in [1, [H, [H, [30J). 

It is clear that, first, if AFi 2 (t) j^ somewhere in two-particle phase space, then AFi 2 (t) ^ almost everywhere, 
since anyway there are many various clusters of collisions resulting in given current states of atoms 1 and 2. 

Second, undoubtedly AFi 2 (t) ^ for post-collision states, since AF\ 2 {t) must compensate and stop nonphysical 
unrestricted propagation of the above mentioned post-collision correlations prescribed by the MC approximation 
F\ 2 {t) = F 12 ■ Hence, undoubtedly AF\ 2 {t) ^ almost everywhere, including the pre-collision states. 

This reasoning shows that appearance of pre-collision correlations along with other ( "non-collision" ) excess corre- 
lations (at states for which AF12 (t) w AF 12 (t) in Eq |3"2"|) is quite inevitable from physical point of view. 

Third, on the other hand, separation of atoms 1 and 2 must decrease a number of the clusters determining C± 2 (t) 
nearly proportionally to visual space angle of one of the atoms from viewpoint of another, ~ ira 2 /4ir\ri 2 \ 2 . The 

matter is that, since F 12 (t) involves fe + 1 (or more) collisions but contains only k integrations, one (or more) of 
k integration velocities is restricted in respect to its direction by a space angle ~ . 
By these reasons, we can propose the following rough fit : 

Ci2(t) = ^—f jj C[ 2 (t) , (33) 

er + 4|r 2 — r\\ A 

where function C 12 (i) keeps non-zero under BGL and smmothly depends on \r\ 2 \ at |r*i2 1 > a • 
Then DVR ([21]) implies (under BGL) that 

' d\r 12 \ f d 3 v 2 C' 12 (t) =™ 2 ^M , (34) 



where = ri 2 /|ri 2 



We see that at any fixed distance |r*i2 1 the excess 
correlations, - particularly, the pre-collision ones, - turn 
to zero under BGL. In this sense, the MC really takes 
place. 

But contribution of these correlations to the "triple col- 
lision integral" ([28]) . to C4 , ... , and thus to the whole 
Eq|3]is determined by region |ri2| ~ a where all excess 
(pre-collision) correlations stay finite under BGL [3l| . 
Hence, in essence MC fails. 

And last remarks. 

(i) Both the left-hand integral in Eq|34land right-hand 
derivative there define some (one and the same) charac- 
teristic length. Of course, it must be nothing but the A 
(for it would be very strange if Cn (t) was indifferent to 
A ). At the same time, obviously, any of constituent parts 
of C 12 (i) , namely, [Z 12 {t) -l\F 1 {t)F 2 (t) and F$(t) 
( k = 1, 2, . . . ), extends up to \r\ 2 \ ~ v$t , with vq be- 
ing characteristic velocity of gas atoms. Hence, indeed 
all these parts are required in order to introduce A in 
place of vot . 

(ii) The reasonings and conclusions of this section (as 
well as that at end of Sec. 5 and in Appendix 2) in no 
way rely on some measure of non-uniformity of gas state. 
Therefore, all they are equally valid, - and thus BE is 
equally invalid, - for arbitrary weakly non-uniform gas 

8. Conclusion. 

Considering spatially inhomogeneous evolution of low- 
density gas, we introduced a non-standard representation 
of its time-dependent distribution functions (DF) in the 
form of their density expansion, and then exploited it to 
derive original exact "dynamical virial relations" (DVR) 
connecting DFs with their density derivatives. Then we 
applied DVR to analysis of behavior of two-particle corre- 
lation function under the Boltzmann-Grad limit (BGL), 
in order to examine validity of the "folklore" opinion that 
evolution of one-particle DF under BGL exactly obeys 
the Boltzmann kinetic equation (BE) while many-particle 
DFs undergo the "molecular chaos propagation" . 

We showed that the corresponding approximate ap- 
proaches to gas kinetics 043 , which seem well grounded 
in the "low-density limit" , at the same time appear non- 
grounded under BGL, since contradict the DVR and thus 
fail, when applied to spatially non-uniform situations (in- 
dependently on degree of the non- uniformity) . 

This fact does not mean, of course, that idea of the 
Boltzmann collision operator (integral) is defective. This 
is excellent concept if one applies it to a separate colli- 
sion. But it by itself is unable to comprise those inter- 
atom statistical correlations what arise from uniqueness 
( "non-ergodicity" (lTj) of histories of collisions in par- 
ticular experiments (see Seel above and comments in 



One of ways to take into account all these correlations 
is the approach to correct solution of the BBGKY equa- 
tions (the "collisional approximation") suggested in [12[ 
(see also [1J, LLJ])- In principle, this approach allows 
to consider a wide variety of phenomena and problems 
( "from molecular Brownian motion to shock waves" [33| ) . 

In [T2I, LLJ] this approach was used to investigate statis- 
tical characteristics of molecular random walks in fluids 
(in particular, the related 1/f noise). The results then 
were confirmed from viewpoint of corresponding exact 
"virial relations" [H, Ell El • Thus, a part of "Augean 
stables" was cleansed: some ancient prejudices, pointed 
out in [llj . like "molecular chaos", were overcome with 
substantial physical profit. 

Therefore, it seems reasonable to try to apply the men- 
tioned approach (strengthened with the DVR) to gas 
(fluid) kinetics and hydrodynamics too. All the more so 
as their difference from Boltzmannian kinetics and clas- 
sical hydrodynamics may be practically important even 
for weakly non-equilibrium (non- uniform) processes [34J . 

Appendix 1. Correlation (cumulant) distribu- 
tion functions and the DVR. 

Let us introduce irreducible many-particle correlation, 
or cumulant, functions (CF) by 



-P12 — F1F2 + C12 , -F123 — F1F2F3+ 
+ C12F3 + C23.F1 + C13F2 + C123 , 



(35) 



and so on. Higher-order CFs can be defined with the 
help of generating functionals: 

00 

S=l 

OO 

= exp[^(o)v] {1 +j2 G Mr/s\} = 

S=l 

OO 

= exp[F 1 (t)i;+J2C s (t)r/sl], (36) 



where, of course, 



F s (t)r = 



Fx... 8 {t)il>{xi)...ip{x s ) 



and so on. In these shortened notations, the DVR (fH))) al- 
together can be accumulated into single generating DVR 



dn 



Sip(x) 



(to be supplemented with equality (J17I) and thus also 
(1201 ). Obviously, this generating DVR is equivalent to 



dn 



Cx. 



. s s+l 



(t) ( s >2), (38) 



s+l 



dn 



= / C 12 (t) (39) 



Generalization of these DVR to dense gases (fluids), - 
when the "initial molecular chaos" (JS]) is too bad choice 
for initial DFs (since they must include effects of atom- 
atom interaction and thus depend on the mean density) , 
- will be considered elsewhere. 



Appendix 2. Many-particle correlations and 
falsity of the BE's "derivations" . 



According to the BBGKY equations (JM) or the DVR 
(TT91 and/or (l3l)|) . the pair pre-collision correlations arise 
in company with various many-particle ones. Their im- 
portance for correct approach to gas (fluid) kinetics was 
demonstrated already in [12| (see also [lj, [lj, IT?} ) • For 
one more demonstration, let us criticize the "derivation" 
of BE (P) suggested in 



Assume, as there, that the irreducible ("pure") three- 
particle correlations (see Appendix 1), - described by CF 
C123 (t) , - can be neglected at a 
BGL. Then the second of Eqs 



n<l and hence under 
1 when written in terms 



L12 G'i2 + L 12 F12 



of CFs, reduces to 

C12 = £12 C12 + L 12 F\ F 2 

( 4 °) 
Solving Eq l40l (with zero initial condition) and insert- 
ing the result, Cn , into the first of EqsJMJ one can 
come (at a 3 n <C 1 ) to the BE Q. 

If this is true derivation of BE, then it should be com- 
patible with the DVR, at least, with Eq]2"Tlfi.e. Eg 155]). 
In fact, however, this is not the case. Indeed, integration 
of EqHH] over x 2 , after multiplying it by n , yields 



li n L Ci2 



L\n\ da 



£'12 ^12 



L\n 



ds 



d_ 



V 



Fi 



where we used also the first of the BBGKY equations 
((Ml) . Combining this equality with the DVR ([21]). we 
come to equality 



[d/dt - Ll][ndF 1 /dn - F x } = 



(41) 



which certainly is wrong. 

Hence, the above BF's "derivation" , based on approx- 
imation (|4"P|) . is erroneous. And, in order to get a correct 
description of iq 's evolution, one should seriously think 
about role of many-particle correlations. 
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